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We study the form factors of vector mesons using a covariant fermion field theory model in 
(3 + 1) dimensions. Performing a light-front calculation in the q + = frame in parallel with a 
manifestly covariant calculation, we note the existence of a nonvanishing zero-mode contribution to 
the light-front current J + and find a way of avoiding the zero-mode in the form factor calculations. 
Upon choosing the light-front gauge = 0) with circular polarization and with spin projection 

h =11= ±, only the helicity zero to zero matrix element of the plus current receives zero- mode 
£Nj \ contributions. Therefore, one can obtain the exact light-front solution of the form factors using only 

the valence contribution if only the helicity components, (h'h) = (++), (H — ), and (+0), are used. 
We also compare our results obtained from the light-front gauge in the light-front helicity basis (i.e. 
h = ±, 0) with those obtained from the non-LF gauge in the instant form linear polarization basis 
5— I 1 (i.e. h = x,y, z) where the zero-mode contributions to the form factors are unavoidable. 
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I. INTRODUCTION 



One of the great challenges in hadronic physics is to calculate the structure of hadrons starting from QCD alone. 
Presently this task is very difficult and one relies on specific models to gain some understanding of hadronic structure 
at low energies and momentum transfer values. A popular model is the constituent quark model (CQM) which in its 
relativized form has met with quite some success. A first test of this model is the comparison of the mass spectra it 
' predicts to the experimental data. Such a test provides some constraints on the wave functions. A more stringent 
test for the wave functions is found when one also calculates the form factors of a hadron. It lies in the nature of the 
£NJ ■ CQM that only valence wave functions are determined easily. However, in a fully covariant calculation of the form 
' factors one needs the full structure of the hadron-quark vertex. 

. It has been known for some time that there are situations where the form factors can be expressed correctly 
O ' as convolutions of the wave functions. Such is the case for certain components of the currents. In particular one 
finds within the formalism of light-front dynamics (LFD) |Q that the so-called plus component of the currents for a 
scalar or pseudoscalar meson can be expressed in terms of the wave functions alone for spacelike momentum transfer. 
The matrix elements obtained this way we call the valence parts. The parts arising from vertices that can not be 
expressed in the wave functions wc call the nonvalence parts. (The plus component of a four-vector being a particular 
combination of its usual components: p + = (p° +p 3 )/v / 2 where the factor V2 is conventional.) 

In the case of vector mesons the situation is more complicated. Till now there have been several recipes for 
the extraction of the invariant form factors from the matrix elements of the currents. It turns out that even when 
one limits oneself to the plus component, these different ways of extracting the form factors do not produce the 
same results j5[. One realizes that since the nine complex matrix elements of the current(J + ), corresponding to the 
possible combinations of the polarizations of the initial and final spin-1 particles, can be expressed in terms of three 
real invariants only, it becomes clear that there must be relations between these matrix elements. This was of course 
known for a long time and many authors have used this knowledge to sort out the invariants from the calculated matrix 
elements. In reference frames where the plus-component of the momentum transfer, q + , vanishes these relations can 
be reduced to just one besides the relations provided by Hcrmiticity, parity and rotation about the z-axis. The latter 
relation is known as the angular condition Q . In general reference frames the situation was not so clear. In a previous 
paper [|| we completely analyzed these conditions for the spin-1 case and found besides the angular condition given 
before another one. There we gave only the formal expressions for these consistency conditions. In the q + = 
frame, however, the additional condition is very simple involving only two helicity amplitudes and doesn't seem to 
provide as strong a constraint as the usual condition since most constituent quark models are expected to satisfy 
it rather easily. Nevertheless, the q + = frame is in principle restricted to the spacelike region of the form factors 
and it may be useful to impose this additional condition in processes involving the timelike region which must be 
analyzed in the q + =/= frame. Thus, it is important to analyze both angular conditions in different frames calculating 
actually the form factors with existing theoretical models. In the present paper, we demonstrate their usefulness for 
thcorctical/phenomenological models for spin-1 objects. In order that the matrix elements satisfy these constraints, 



the current operator must transform properly and the state vectors must be eigenstates of total spin. If the models 
do not have these properties, the angular conditions will not be met. 

In this work, we use a simple but exactly solvable model for the spin-1 (e.g. p) meson and separate the valence 
and nonvalence contributions to the three physical form factors to investigate the degree of violations in the two 
angular conditions for each contribution in different frames. Although the quantitative results that we find in this 
model may differ in other models depending on the details of the dynamics in each model, the basic structure of 
model calculations is common and we expect the essential findings from this model calculation may apply to realistic 
models. 

In particular, we compared two different types of polarization vectors, the one obtained from the LF gauge = 
0), which is usually used in the LF CQM analysis, and the other obtained from the instant form (IF) polarization, 
which is not associated with the LF gauge, i.e. e+ ^ 0, but used in some recent papers (?] |]]. In both cases, there 
is a zero-mode contribution, i.e. a contribution from the nonvalence part that remains finite for q + — ► 0, even if 
the plus-component of the currents is used. Specifically, there is a zero mode contribution in the LF helicity case 
(h = +,—,0) to the (h',h) = (0,0) amplitude, where h and h! are the initial and final helicities, respectively, but 
there is no zero-mode for other helicity combinations such as (+,+), (+, — ) and (+,0). On the other hand, in the 
instant form case, only (yy) is immune to the zero- mode but others such as (xx), izz) and (zx) do receive zero- mode 
contributions. Of course, the two results are exactly the same if one properly includes the zero-mode contribution. 

Now turning to the angular conditions, there are several different prescriptions in choosing the matrix elements 
to extract the three physical form factors. We compare three different types of helicity combinations, GK 0], CCKP 
||, and BH using both LF and instant form helicity bases in a reference frame where q + = 0. One of our very 
interesting findings of the analysis in the LF helicity basis is that the prescription using the plus-component of the 
current but not involving the (h',h) = (0,0) helicity amplitude in the LF gauge is preferred for model calculation. 
Especially, the GK prescription uses only (hi , h) = (+, +), (+, — ) and (+, 0) but not the pure (0, 0) component and thus 
achieves the goal of not involving the zero-modes. On the other hand, the longitudinal (0,0) component is the most 
dominant contribution in the high momentum transfer region and thus it may be better to use the BH prescription, 
involving the (0, 0), (+, 0) and (+, — ) amplitudes only, in the high momentum perturbative QCD analysis. The CCKP 
prescription, however, involves all helicity states, i.e. (+,+), (0, 0), (+, 0) and (+, — ) and one needs a quantitative 
analysis of the angular conditions to pin down the momentum transfer region for the validity of this prescription. 
Our quantitative analyses indeed verify that the GK prescription is remarkably free from the zero-mode contribution 
but others are not. In the recent work by Melikhov and Simula we see that the result using the GK prescription 
is not in complete agreement with their covariant model calculation, which is due to the dependence on a light-like 
four-vector called w(w 2 = 0) in their formulation that necessarily involves unphysical form factors. The covariant 
formulation presented in our work should be intrinsically distinguished from theirs Q , since our formulation involves 
neither u> nor any unphysical form factor. 

If we use the instant form basis, however, all three prescriptions receive the zero-mode contribution. Since the 
instant form helicity is not obtained from the LF gauge, i.e. e+ ^ 0, even the GK prescription gets the zero-mode 
contribution especially for the magnetic form factor as we shall show in this work. Thus, the instant form basis used 
in the LF formulation seems quite dangerous, because it can lead to a wrong interpretation of the physics involved in 
LF dynamical models. Our solvable model calculation clearly indicates that one can avoid the zero- mode contribution 
if the LF basis when the LF gauge is used without using the longitudinal to longitudinal helicity amplitude. 

This paper is organized as follows. In Section H, we summarize the angular conditions for spin-I systems using 
the LF helicity basis and the kinematics for the reference frames Drell-Yan-West (DWY), Breit (BRT), and target- 
rest frame (TRF) used in this work. The three prescriptions (GK, CCKP, BH) used in extracting the physical form 
factors are also briefly discussed in that Section. In Section |TJ, we present our covariant model calculations of physical 
quantities such as the three electromagnetic form factors and the decay constant of the spin-1 meson system using 
both the manifestly covariant Feynman method and the LF technique. In the q + = frame, we separate the full 
amplitudes into the valence contribution and the zero-mode contribution to show explicitly that only the helicity 
zero to zero amplitude is contaminated by the zero-mode. In Section [iv| , we present the numerical results for the 
form factors and the angular conditions and analyze the dependences on the prescriptions, reference frames, and 
helicity bases. The taxonomical decompositions of the full results into valence and non- valence contributions are used 
wherever possible to make a quantitative comparison of these dependences. Conclusions follow in Sectio n M . The 
details of the instant form analysis and a derivation of the zero-mode are summarized in the Appendices |A] and |^, 
respectively. 
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II. SPIN-1 FORM FACTORS IN LIGHT-FRONT HELICITY BASIS 



The Lorentz-invariant electromagnetic form factors -Fx, F 2 , and F3 for a spin-1 particle of mass m are defined flCf l 
by the matrix elements of the currents J M between the initial \p, h) and the hnal \p' , h') eigenstates of the momentum 
p and the helicity h as follows: 

Gh'h = (P', h'\J"\p, h) = -el ■ e h (p + pYFi(Q 2 ) + (e£q- e* h , - e* h f q ■ e h )F 2 (Q 2 ) + ^^^ (P + p'TMQ 2 ), 

(1) 

where Q 2 = —q 2 , q — p' — p and thi^h') is the polarization vector of the initial(final) meson. 

The physical form factors, charge, magnetic, and quadrupole, are related in a well-known way to the form factors 
Fi, viz, 

Gc = (l + \r,)Fx + \vF 2 + + V )F 3 
Gm = — F2 

G Q =F 1 + F 2 + (1 + V )F 3 , (2) 

where 77 = Q 2 /(4m 2 ). 

Using the convention e M = e~ , e 1 , e 2 ), the general form of the LF polarization vectors is given by 



sMp + ,p\p 2 ,o) = fe 2 ^,^)- (3) 

e^ F (p+,p 1 ,p 2 ;-) J ! / pi_ j_ _i> 

P +> 72 > 72, 





G+o 


GX 




r+ 


Gq 


G± + 


G_ 


G± 



Here p r (p l ) = T(Px ± iPy)/V^- Using Eqs. (|T|) and (||), we obtain the matrix elements: 

Gt'h = 

\G_ 

(aiFt + a 3 F 3 ci.Fi + c 2 F 2 + c 3 F 3 e 3 F 3 

biFi + b 2 F 2 + 63F3 diFi + d 2 F 2 + d 3 F 3 -(b t Fi + b 2 F 2 + b s F 3 )* ). (4) 

e 3 F 3 -(ci-Fi + c 2 F 2 + C3F3)* 01F1 + a 3 F 3 

Since we are working only with the plus component of the current, we shall use the following short-hand notations 

a — ( - T ++ — ^ , ^6 — ^0+ ~~ — ^O— ' 



G c — Ga 



f~l+* 



G d 



G e = G 



G 



+* 



(5) 



The invariant form factors can be extracted in a straightforward way. The simplest procedure is to solve hrst for 
F3 from G e - Next F\ is obtained from G a and F 3 . Then there are three options for obtaining F 2 from G&, G c , and 
Gd- These solutions are denoted by Fj* > -^21 an( i ^2 > respectively. The full solutions are then 
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ai 
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G, 
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G c 



Ga 



a 3 b x - aib 3 

a x e 3 
«3Ci - QlC 3 

a\e 3 
a 3 d\ - aid 3 

a\e 3 



G, 



G, 



G, 



(6) 
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This procedure makes it clear that the covariant form factors of a spin-1 hadron in Eq. ([j]) can be determined using 
only the plus component of the currents, G^, h (0) = (P 1 , h'\ J M |P, h) 7 in any chosen Lorentz frame. The nine elements 
of the current operator G^, h (0) must be constrained by the invariance under (i) time-reversal, (ii) rotation about z 
and (iii) reflection in the plane perpendicular to z, and rotational covariance, i.e. invariance under the rotations about 
an axis perpendicular to z. So two additional constraints on the current operator are required. These consistency 
conditions are the angular conditions, which we define as |^| 

A bc = F| - F 2 C , A bd = F|-F 2 d . (7) 

We know that the form of these conditions depends on the reference frame m . In this work we consider three different 



frames, which we define in the following subsection [I A. Especially, in the frames where q + = 0, our angular condition 
Abd is equivalent to the usual angular condition relatin g th e four helicity amplitudes discussed in the literature 0] 



modulo an overall factor as we discuss in the subsection II B 



A. Kinematics 

Our conventions for the momenta of the initial and final state mesons in the three different reference frames, 
(DYW), (BRT), and (TRF) are given below. We use the notation p M = (p + ,p~ ,p x ,p y ) = (p + ,p~,p±) and the metric 
convention p ■ q — p + q~ + p~q + — p±_ ■ qx- 



p=(p+,m 2 /(2p+),0 1 0) 

p' = (P + ,(Q 2 + m 2 )/(2p+),Qcos0,Qsm<f>) (8) 




TRF 



DYW 



BRT 

/2m/3 — Q cos 9 2mf3 + Q cos 9 Q sin9 cos <fi Qsin0sin< 
P= \ V2 ' V2 ' 2 ' 2 
. / 2m(3 + Q cos 9 2m(3 — Qcos9 Q sin 9 cos <f> Qsin6>sin0\ 
P= { ' 2~T2 ' 2 ' 2 J' (10) 

(11) 

(mm 

, fm + k + PQcosd m + k- f3Qcos9 . . . \ 

p = y , ,PQsm0cos(p,pQsm9sm(pj. (12) 

In the literature usually the reference frames used are limited to ones where q + = (q 2 = 2q + q~ — q\ < 0). One of 
such reference frames is the special Brcit frame used in Refs. |2]-|^, |ll| , |l2| , where q + — 0, q y = 0, q x = Q, and p± = —p'j_ 
i.e. 6 = f,</> = in Eq.®; 

Special Breit 

f = (0, 0, Q, 0), p» = (m^/TT^/V^ m^TT^/V^ -0/2, 0), p'» = (m^TT^/V^ m^TT^/V^ 0/2, 0), (13) 

where 77 = Q 2 /Am 2 is a kinematic factor. The corresponding polarization vectors are obtained by substituting these 
four vectors in Eq. (^|) and the transverse(/i = ±) and longitudinal (h — 0) polarization vectors in this special Breit 
frame are given by 
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eV ' ±) ^( '^ 1 ' ±l )'^' 0) = ^' 2;, 



74 



(14) 



B. Angular condition in q + = frame and prescriptions of choosing helicity amplitudes 



In the q + = frame, one can reduce the independent matrix elements of the current down to four, e.g. 
G+ + ,G+_,G^ and G^ using the front-form helicity basis [p|-[5|,|l2"t and the rotational covariance requires one addi- 
tional constraint on the current operator. This is what these authors call the angular condition A(Q 2 ) and can be 
obtained from the explicit representations of the helicity amplitudes in terms of the physical form factors. Using the 
relation between the covariant form factors F{ and the current matrix elements given by Eq. (|l]), one can obtain the 
following helicity amplitudes in the q + = frame 



G+ + = 2 P +(F 1 + 77F3), G+ = p+^{2F 1 +F 2 + 2 V F 3 ), 
G+_ = -2p+ V F 3 , G+ = 2p+{(l - 2 V )F 1 - 2r,F 2 - 2r? 2 F 3 }. 



(15) 



Thus, the usual angular condition relating the four helicity amplitudes is given by 



A(Q 2 ) = (1 + 2 V )G+ + + G+_ - y^G+ 



Gnn 



0. 



(16) 



where we note an overall factor difference between A(Q 2 ) and A^Q 2 ), i.e. A = c^A^ (See Section [VB5 for the 
discussion of the factor d.2-). 

In a practical computation, instead of calculating the Lorentz-invariant form factors Fi(Q 2 ), the physical charge 
(Gc), magnetic (Ga/), and quadrupole (Gq) form factors are often used]]. However, the relations between the physical 
invariant form factors and the matrix elements G~^, h || are not unique. Only if the matrix elements fulfill the angular 
condition Eq. ( |l6| ) the extracted form factors would not depend on the choice rn.ade._So one may choose which matrix 
elements to use to extract the form factors. Perhaps the most popular choices are 



GK 

c — 


1 

2p+ 


[(3-27?) 
3 


G+ + + - 


GK 
M — 


2 
2p+ 


G + ++ - : 
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(17) 



G CCKP 



pCCKP 



Gq 



1 
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2p+(l + »7) 
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3 y^j ' 


6 


2(l-r?)„ + " 

V2H G+0 _ 






Gio fa + 2) 
2t? 
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(18) 



and 



x In Refs. |^J^, the form factors Go, Gi, and G2 are used and the two definitions are related by Gc ~ Go/2p + , Gm = Gi/2p~ i 
and (t7v / 8/3)Gq = G 2 /2p+ 
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■^BH _ 1 

c 2p+{l + 2 V ) 

r" BH = - 

2p+(l + 2 V ) 



1 



2p+(l + 27?) 



(3-2t?) + 

^ Crnn 



_ fi-\ 



1677 G+ 2(277-1) 



3 ^ 
(277-I) 



2r] G+0 



9 G +0 ^+ (1 + 7/) ^ 



G 



(19) 



The relation between F's and G's given by Eq. (g) holds for any prescription given above. 

It is interesting to note that while Grach and Kondratyuk in Q regarded Gq as the "worst" element and took 
care not to use it writing the relations Eq. (0), Brodsky and Hiller Q included Gq instead of G+ + expecting the 
helicity zero to zero component of the current matrix element to be the dominant one in the perturbative QCD regime. 
Chung et al. j^] used all four independent helicity components of the current matrix elements. On the other hand, 
in the instant form basis used by some authors |7|fj|,|l3[] , the independent matrix elements of the current operator are 
G+j., Gy V , G + z and G+, and the angular condition becomes A(Q 2 ) = Gy y — G+ z . In appendix A, we show the relevant 
expressions for the form factors in the instant form spin basis Although the authors in Refs. §,|,|l3| argued 

that this basis is completely equivalent to the LF helicity basis, the relation between them, which will be discussed 
later, is not trivial. 



III. CALCULATION IN A SOLVABLE COVARIANT MODEL 



The solvable model, based on the covariant Bcthc-Salpeter(BS) model of (3 + l)-dimcnsional fermion field theory, 
enables us to derive the form factors of a spin-1 particle exactly. 




FIG. 1. The covariant triangle diagram (a) is represented as the sum of a LF valence diagram (b) denned in the region 
< k + < p + and the nonvalence diagram (c) defined in p + < k + < p' + . 5 — q + /p + — p' + /p + — 1. The white and black blobs 
at the meson-quark vertices in (b) and (c) represent the LF wave-function and non-wave-function vertices, respectively. The 
small circles in (b) and (c) represent the (on-shell) mass pole of the quark propagator determined from the ^"-integration. 



The covariant diagram shown in Fig. |l|(a) is in general equivalent to the sum of the LF valence diagram (b) and the 
nonvalence diagram (c), where S = q + /p + = p' + /p + — 1. The matrix element G^, h (0) of the electromagnetic (EM) 
current of a spin-1 particle with equal mass constituents (m q — rriq) obtained from the covariant diagram of Fig. |l](a) 
is given by 

C» (0) =M a 2 ( — S A (k-p)SZ, h S A {k-p') 

h ' hK 1 cy J (2n)± [(k - p) 2 - m 2 q + ie][k 2 - m 2 q + ie][(k - p>) 2 - m 2 q + ieY K ' 

where N c is the number of colors and g, modulo the charge factor e q , is the normalization constant, which can be 
fixed by requiring the charge form factor to be unity at zero momentum transfer. Sy h is the trace term of the quark 
propagators. To regularize the covariant fermion triangle-loop in (3 + 1) dimension, we replace the point photon- vertex 
7 M by a non-local (smeared) photon-vertex Sa(p)7 m Sa(p'), where Sa(p) = A 2 /(p 2 — A 2 + is) and A plays the role of 
a momentum cut-off similar to the Pauli-Villars regularization [jjj . 

When we do the Cauchy integration over fc~ to obtain the LF time-ordered diagrams, we want to avoid the 



complexity of treating double k -poles, so we decompose the product of five energy denominators in Eq. (20) into a 
sum of terms with three energy denominators only: 



1 1 1 
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\D'a 


D'J 



(21) 



G 



where 



D A =(k- pY - A z + ie, D = (k- pY -m\ + ie, D k = k l - m\ + ie, (22) 
and D' Q[A] = D 0[A] (p -> p'). 

Our treatment of Sa as the non-local smearing photon- vertex remedies |jj the conceptual difficulty associated with 
the asymmetry appearing if the fermion-loop were regulated by smearing the qq bound-state vertex. As will be 
discussed later, the two methods lead to different results for the calculation of the decay constant even though they 
give the same result for the form factors. 

The vector meson decay constant /y in this covariant model with the nonlocal gauge boson vertex S\(k)^f fl SA(k — p) 
is defined by 

A" = (0\q^q\p; 1J 3 ) = tV2f v me^(J 3 ), (23) 

where 

AM = N ,4 f _d*k Tt[^- p'+m q )r(fi + mq )} 

cJ J (2tt) 4 [k 2 - m 2 + ie][(k - p) 2 - m 2 + ie][k 2 - A 2 + ie][(k - p) 2 - A 2 + ie]' 1 ' 



A. Manifestly covariant calculation 

In the manifestly covariant calculation, we obtain the form factors Fi(i = 1,2,3) using dimensional regularization. 
Although the splitting procedure Eq. ( pl| ) may not be neccessary in the covariant calculation, it seems more effective 
in practical computation. Here we describe some essential steps for the derivation of the covariant form factors: We 
(i) reduce the five propagators into the sum of three propagators using Eq. (^l]), (ii) use the Feynman parametrization 
for the three propagators, e.g., 



= 2 / dx 

'h. D D Q Jo Jn 



dy 



[D k + (D Q - D k )x + (D' Q - D k )yf 



(25) 



and (Hi) make a Wick rotation of Eq. ( p0| ) in D-dimension to regularize the integral, since otherwise one encounters 
missing the logarithmic divergent terms in Eq. (po|). Following the above procedures (i) -(Hi) we finally obtain the 
covariant form factors as follows: 



F!(Q 2 ) = 



F 2 (Q 2 ) 



F 3 (Q 2 ) 



N c g 2 A 4 



1 pl-x 

dx / dy{ (2 
i) Jo 



y)\n 



<~<2 f^2 
^fcAO^fcOA 

^kAA^kOO 



8tt 2 (A 2 „.„ 

-(x + y)(x + y — l) 2 m 2 + (2 - x — y)xy Q 2 - (2 - x - y)m 2 



c- 



N c g 2 A 4 



dx 



dyU2 



y)ln 



r 2 r 



fcAO^fcOA 



8^ 2 (A 2 ~m 2 ) 2 7 Jo 
(x + y)[(x + y) 2 - l]m 2 + (x + y)xy Q 2 - (2 + x 
N c g 2 A 4 



i~i2 r~i2 
°feAA°fcOO 



y)m 2 q 



8tt 2 (A 2 -m 2 ) 2 J 



dx 



dy 8xy(x + y — l)m C 



(26) 



where 



(i/c, 



^feAA = + y)(l -x- y)m 2 -xyQ 2 - (x + y)A 2 - (1 



y)r 



CfeAo = (x + y)(l -x- y)m 2 - xy < 

CfeOA = CfeAo( x 2/)' 

(x + y)(l — x — y)m 2 — xy < 



(xA 2 + ym 2 ) -(\-x- y)m 2 , 



r 2 

^kOO 



and C 2 

the dimensional regularization. 



fcAA 



1/Cf 



feAO 



(27) 



feOA 



1/Cf 00 ). Note that the logarithmic terms in F± and F2 are obtained from 
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Following a similar procedure for the form factor calculation, the covariant result for the decay constant is obtained 

as 

/-v = ^ 4 f d J H + ^ ^[^-^-(l-^y-^-d-.K] 

4\/27r 2 m(A 2 - m2)2 7 \ L « [M 2 — A 2 ] [M 2 — m 2 ] 

— [M 2 - A 2 ] ln[-M 2 + A 2 ] - [M 2 - m 2 } ln[-M 2 + m 2 ] 

+ [M 2 - xm\ - (1 - x)A 2 ] ln[-M 2 + xm\ + (1 - x)A 2 ] 

+[M 2 - 2A 2 - (1 - x)m 2 ] ln[-M 2 + a;A 2 + (l-a;)m 2 ]|, (28) 

where M 2 = x{l — x)m 2 . 



B. Light- front calculation 



We shall use only the plus-component of the current matrix element Gtrt, in the calculation of the form factors. In 
principle, one can directly calculate the trace term S£, h with k~ = fc~ ole , which depends on the integration region of 
k + . However, for the purpose of a clear understanding of the physics implied in LF dynamics, we instead separate 
S^, h into the on-mass shell propagating part and the (off-mass shell) instantaneous one using the following identity 



+ m q )+~ 1 + {p -p on ), 



(29) 



where the subscript (on) denotes the on-mass shell (p — m 2 ) quark propagator, i.e. p = p on = (to 2 +p 2 L )/2p 



Then the trace term S~^, h of the quark propagators in Eq. (|20| ) is given by 

St h (P =k-p,P> = k-p>,k) = (5+ Jon + (S&fcW, 

where 

(S+ h ) on = TrlfiAKn + ™ 9 )7 + ( An + m q ) + m,)], 

(L • Cfc') (-Pon ' Eft) + {Pon ' 4')( fc on ' + {f* h , ' £h){m - k 



(30) 



4 P+ 
+ 4P'+ 
+ 4 fc+ 
"4^ 



(fc on ' 4') {Pon ■ £h) + {Pon ' £*h>){ k on ' £/t) + (4< ' e/ l )(^ 2 - fc on ■ Pon) 

(P on • e£,)(P on • e h ) - (Pon • e* h ,){PU ■ e h ) - (e* h , ■ e h ){m 2 - P on ■ PL) 



{P on ■ e h ){k on ■ P on - m 2 ) - (P^ n ■ e h ){k on ■ Pon - ml) + (fc on • e h )(p on ■ P' on - m 2 ) 



4 4 



{Pon ' e h'){kon ' Pon ~ m„) — (P on ' e h ,){k on ' P +on ™ ? ) + (fcon ' £h')(Pon ' P on — m o) 



(31) 



and 



(5+ Ji„st. = (ft" - fcon)TrK,(P on + m q ) 1 +{ P on + m,) rf, 7 + 
= 8(fc- - fe- ) 



^P Q +(e h • P on ) + e+P '+(e^, • P on ) - P +P +(e^, • e h ) + e£+e+(m 2 - P' on ■ P on ) 



(32) 



As we shall show below, the LF valence contribution comes exclusively from the on-mass shell propagating part, 
Eq. (|3l"|), and the zero-mode (if it exists) from the instantaneous part, Eq. (|32j). 

Using the special Breit frame (See Eq.(|l3|).) with the LF gauge, we obtain for the trace terms {S^, h ) on and {S^, h )i ns t. 
given by Eqs. (|3ll) and fl3^ ) the expressions 



(^++)° n 



4p 4 

a; 



to 2 + (2x 2 - 2x + 1) I fci - — Q 2 + ix{k± x <fj_) ■ z 



8 



and 



(»?+ — )on 
("^+o)on 
(S'oo)on 



8(1 - x)p~ l 



Q 

Ap+ 



p + (2(k x - iky) - xQ)(2x - 1) 



x[ (l-s)(m a + M 2 ) + -Q 2 



(1 - x){m 2 + M 2 ) + jQ 2 +k X -q± 



- x(k± ■ q±f 





inst. 


— )inst. — 


o, 






(£f0 


inst. 


= 8 <» + >V- 


- Kn)( x - 1) 






(^00 


inst. 




r x 2 


2 + m 2 + k\ 





(33) 



(34) 



where x = k + /p + , and Mg = (m 2 + kj_)/[x(l — x)]. We note that the terms proportional to an odd power of kj_ do 
not contribute to the integral. 

By doing the integration over k~ in Eq. J2G|), one finds the two LF time-ordered contributions to the residue 
calculations corresponding to two poles in k~, the one coming from the interval (I) < k + < p + (See Fig. |I](b).), the 
"valence diagram", and the other from (II) p + < k + < p' + (See Fig. |l](c).), the "nonvalence diagram" or "Z" graph. 
These diagrams are expressed in terms of energy denominators. 



1. Valence contribution 

In the region < k + < p + , the pole k~ — k~ n = (m 2 + k\ — ie)/2k + (i.e the spectator quark), is located in the 
lower half of the complex fc~-plane. Thus, the Cauchy integration formula for the ^"-integral in Eq. ( po|) gives in this 
region for the plus current, G^, h (0), 

r +val _ N c f 1 dx f ,2? 9^ q+ val 9_^_ 

h ' h 2(27r)3i x(l~x) 4 J ^{m?-M 2 ){mZ-M\) h ' h (to 2 - M l2 ){m 2 - M%Y 1 ' 

where 



m 2 + (k± - xp±) 2 to 2 + (k ± - xp ± ) 2 



1 -x 



M 2 A = m i + (^-^) 2 + ^ 2 + (^-^±)\ (36) 

X 1 — X 

are the invariant masses of the initial meson state. The invariant masses of the final state, i.e. A4'q and Ai'j 2 in 



Eq. (pa), can be obtained by replacing p^ — > — p± in Eq. (36). As one can easily see from Eqs. (J33|) and (134) , only 



the on-mass shell quark propagator part contributes to the valence diagram, i.e. S^™ — (S£, h ) on . Note, however, 
that this relation does not hold in general for other components of the currents, e.g. S^,™ ^ (S^, h )on- One of the 
distinguished features of the LF plus current matrix element given by Eq. (|3^) is that the physical interpretation is 
manifest in terms of the LF wave function, i.e. a convolution of the initial and the final state LF wave functions, 
which is not possible for the covariant calculation. 



2. Zero-mode contribution 

In the region p + < fc+ < p' + , the poles are at k~ = jf~ + [m 2 + (k± —p'±) 2 — ie]/2(k + — p' + ) (from the struck quark 

propagator) and = j/~ + [A 2 + (k± — p' j_) 2 — ie]/2(k + — p l+ ) (from the smeared quark-photon vertex S\(k —p')), 
and are located in the upper half of the complex fc _ -plane. 

Since the integration range of the nonvalence region, p + < k + < p' + (= p + + q + ), shrinks to zero in the q + — ► 
limit, the nonvalence contribution is sometimes mistakenly thought to be always vanishing for q + — > 0. However, in 
reality it may not vanish but give a finite contribution, 
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fP + +q + pl+8 

lim / dk + (---) = lim / dx(- • •) ^0 



(37) 



Then it is called the "zero-mode" [jl4j-|l9j in the q + = frame. The nonvanishing zero-mode contribution occurs only 
if the integrand (• • •) in Eq. ( |37| ) behaves ~ (i.e. (1 — Note that there is no zero-mode contribution either in 

the case the integrand behaves like k~(k + —p + ) n (n > 1) or is fc~-independent. 

For the plus current, the zero-mode contribution comes from the spin structure of the fermion propagator, specifically 
only from the instantaneous part given by Eq. (|34| ) and neither from the on-mass shell propagating part nor the energy 
denomimator. Thus, without detailed knowledge of the energy denominator, it is easy to find from Eqs. ([54|) and ( [37] ) 
that only the helicity zero to zero component gives a nonvanishing zero-mode contribution: 



lim (5+ 



00/inst. 



8(p+) : 



k (m 



74), 



(38) 



where k~ ~ 1/(1 — x) — ► oo as x — ► 1. In other words, while the integration region shrinks to zero, the integrand for 
the helicity zero to zero component goes to infinity leading to a finite zero mode contribution. 

As we said before, we avoided the complexity of the Cauchy integration over double /c~-poles by decomposing the 
product of five energy denominators in Eq. ( |20| ) into a sum of terms with three energy denominators. In this way, we 
perform the Cauchy integration of Gq over the single fc~-pole, either D' A or D' , instead of double fc~-poles. 

For example, the l/{DkD\D' A ) term in Eq. (pl|)_combined with the pole position fc~ = p'~ + [A 2 + (fci— p' x) 2 ]/2{k + — 
p' + ) appearing in (S' c J )i ns t. gives (See Appendix [b| for the detailed derivation.) 



d A k 



{2ixY D k D A D' A 2(2tt)V 



In 



d 2 k. 



(k ± -p' ± ) 2 +A 2 



[(fcx - Px) 2 + A 2 ] - [(fc x - p'x) 2 + A 2 ] 



(39) 



Similarly, one can obtain a nonvanishing zero-mode contributions for the other energy denominator terms given by 
Eq. (Elf). Explicitly, the zero- mode contribution from Sq in this special Breit frame is given by 



N c g 2 A A 



2p+(2tt) 3 (A 2 -m 2 ) 2 



,2t 8 (P + )\ 2 

d k± 5— (m 



kl 



2 /4) 



In 



(fc^-p'^) 2 +A 2 

(fe_L-pjJ 2 +A 2 



hi 



l(k± - pi) 2 + A 2 ] - [(k± - p\Y + A 2 ] [(k± - p±) 2 + A 2 ] - [(k x - p'j 



+ {m q ^A)\. (40) 



The angular condition A(Q ) given by Eq. (nm is satisfied only if the zero- mode contribution for Cgg m ' in Eq. (W 



is included, i.e. G> 



G, 



+ val 
00 



^OO 



A similar analysis has been made by de Melo et al. M, where the authors found the zero- mode contribution using 
the instant form basis ||],|],[l3| instead of the LF helicity basis [p[-p| JTo| , ^2] | for the polarization vectors of a spin-1 
particle. In principle, the LF helicity basis can be related to the instant form spin basis by some transformation. 
Interestingly, however, we find that since the authors in Ref. used the non-LF gauge(i.e. e+ 7^ 0) polarization 
vectors, the three polarization components, i.e. G^ X ,G^ Z and Gj x , receive zero-mode contributions as we explicitly 
show in appendix A. In other words, using the instant form basis with a non-LF gauge used in one cannot avoid 
the zero-mode contribution to the form factors of a spin-1 particle no matter what prescription is used. 

We use the results of our numerical calculations to compare the form factors obtained in the LF helicity basis (in 
LF gauge) with those obtained in the instant form linear polarization basis (in non-LF gauge) as well as the covariant 
ones. 

In the LF calculation of the vector meson decay constant, the plus current with the longitudinal(/i = 0) polarization 
vector is usually used. In the special Breit frame (See Eqs. JT^) and ([l4|).), we thus obtain 



LF 
V 



N c gA 4 
4V2~7r 3 r 



dx 



x 3 (l — x) 3 

2x(l — x)m 2 



d 2 k. 



x(l - x)(p + Y + m\ + kj_ - kx ■ px 
) - A 2 - 2(fcx - xpx) 2 



[to 2 - M 2 J[to 2 - M 2 A ][m 2 - Mf )m {m q <-> A)][m 2 - M 2 A (rn 9 «-> A)] ' 



(41) 



Our LF calculation of the decay constant in Eq. (41) is exactly the same as the covariant result in Eq. (|28|). We also 
note that there is no zero mode contribution to fy in our model calculation. This can be easily seen from the trace 
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calculation, because Tr[y($— fl+ m q )j + (ft + m q )] =4{(e • k)(2k + — p + ) + e + (m 2 — k 2 + k-p)} and the fc~-terms cancel 
each other, so only the good component is left in the numerator. It is interesting to note that while our calculation 
of the decay constant with a non-local (but symmetric) gauge boson vertex is immune to the zero-mode, the same 
calculation by Jaus |2p| is not, where the author used a local gauge boson vertex and an asymmetric smearing meson 
vertex. 



IV. NUMERICAL RESULTS 



In this section, we present the numerical results for the form factors and angular conditions and analyze the 
dependences on prescriptions, helicity bases and reference frames. However, we do not aim at finding the best-fit 
parameters to describe the experimental data of the p meson properties. Rather, we simply take the parameters used 
by others [g| with which we were able to reproduce the results in that particular work. Nevertheless, as we mentioned 
earlier, our model calculations have a generic structure and the essential findings from our calculations may apply to 
the more realistic models, although the quantitative results would differ in other models depending on the details of 
the dynamics in each model. 

In our numerical calculations, we thus use m = 0.77 GeV, m q = 0.43 GeV, and A = 1.8 GeV || and make the 
taxonomical decompositions of the full results into the valence and nonvalence contributions to facilitate a quantitative 
comparison of the various dependences such as the prescriptions (GK,CCKP,BH), the helicity bases (LF,IF) and the 



reference frames (DYW,BRT, TRF). We first present the de 
q + = frame (See subsection IV A[ ). Then, in subsection 
the LF helicity basis. 



jcnde nces on the prescriptions and the helicity bases in the 



IV B, we present the frame dependences using exclusively 



A. Dependences on the helicity bases and the prescriptions 

In Fig. ^, we show the charge form factor \Gc(Q 2 )\ obtained from the light-front (left) and the instant-form (right) 
spin bases. The full solutions (thick sol id li ne) are obtained from three different prescriptions |2[jJ] given by Eqs. (|l7j)- 
( |l9| ) for the light-front basis and Eqs. flAl| )-([A"3]) for the instant-form basis, respectively, and they all turn out to give 
exactly the same result as the covariant one as they should be. The slope of the full solution gives the charge radius 
of the bound-state as defined in Eq.(H) and we obtained (r 2 -,) = 7.63 GeV~ 2 with the param eter set we used. More 



detailed discussions on the charge, magnetic and quadrupole radii can be found in subsection IV B 6. In the q + = 
frame, the full solutions can be decomposed into the valence contribution and the zero-mode contribution since the 
nonvalence diagram reduces in the limit q + — > to the zero mode. To estimate it, we plot the valence contribution 
for each prescription, i.e. the dotted line for GK |Q, long-dashed line for CCKP [Bj, and dot-dashed line for BH [Q, 
respectively. The normalization constant g is fixed by requiring the full solution to be normalized to Gp(0) = 1. As 
one can see in Figj2j the two results for the valence contributions obtained from the light-front and the instant-form 
bases exactly coincide with each other. However, only the GK prescription is immune to the zero-mode contribution 
for both helicity bases. The dotted curve cannot be seen because it is on top of the solid curve. Other prescriptions, 
CCKP and BH, receive large amounts of zero-mode contributions (i.e. the difference between the full solution and 
the valence one). As we discussed earlier, the GK prescription does not involve the Gq component which is the only 



source of the zero-mode for the light-front helicity basis and the zero-modes from the G xx and r]Gf z terms in Eq. (Al) 



for the instant-form basis cancel each other(See Eq. (|A6|).). We also show the angular condition(small squares) given 



by Eq. (16) without including the zero-mode contributions. If we include the zero-mode contributions, then it is of 
course exactly zero. 

The situation is rather different for the calculation of the magnetic form factor Gm as shown in Fig. ^. For the 
full solution, the two (LF and IF helicity bases) results are again exactly the same as they should be. The magnetic 
moment (i n units of e/2m) and its radius given by Eq.(44) are obtained as /xi = 2.1 and (r 2 M ) = 9.73 GeV" 2 (See also 



subsection [VB6] .), respectively. However, the valence (or for that matter the zero-mode) contributions to the full 



solution are quite different depending on the helicity bases. For the light-front helicity basis, the GK prescription is 
again immune to the zero-mode and the dotted curve is exactly on top of the solid curve. Also, the other prescriptions, 
CCKP and BH, receive large amounts of the zero-mode contributions as in the case of the Gc calculation. However, 
for the instant form spin basis used in f|, not only the CCKP and BH prescriptions but also the GK prescription 



are affected by the zero- mode, because the zero- mode terms — G+ and G* x / \/rj in Eq. (Al) do not cancel each other 
(See Eq. (|A6|).) but rather add up. 

We show in Fig. [| the quadrupole form factor Gq{Q 2 ) obtained from the light-front (left) and the instant-form 
(right) helicity bases. The quadrupole moments (in units of e/m 2 ) and the corresponding radius given by Eq. (fl4|) are 
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obtained as Qi = 0.91 and 



Qi 



12.6 GeV 2 (See also subsection IVB 6.), respectively. As in the case of Gc(Q 2 ), 



the two (LF and IF hclicity bases) results coincide and the dotted curves are exactly on top of the solid curves because 
of the absence of the zero-mode in the GK prescription. 

The decay constant (See Eqs.(|2"8|) and pi]).) using the same parameters yields the result f v = 133.7 MeV, while 
the experimental data f p o = 152.8 ± 3.6 MeV and f p ± — 147.3 ± 0.7 MeV are obtained from the width T(p — > e + e~) 
and the branching ratio Br(r — > pv T ) = (25.02 ± 0.16)% Q, respectively. 



B. Light- front valence parts 



We checked that in all reference frames the sum of the valence and nonvalence contributions to the form factors 
is equal to the covariant result. Therefore henceforth we plot the valence contributions only. The valence parts will 
in general depend on the polar angle 9 in BRT and TRF, but are independent of the azimuthal angle ip in all three 
reference frames(DYW,BRT,TRF). We used the latter property as a check of the accuracy of our codes. 

Using Eq. (0) for the matrix elements and the kinematics specified in Eqs. ([Io]) and ( |l2|) for BRT and TRF, 
respectively, one finds that the coefficients hi and Ci, i = 1,2, 3, vanish for 9 — 0. Therefore we illustrate the angular 
dependence of the valence parts in BRT and TRF by giving them for the small but nonvanishing value 9 — 7r/20. On 
the other hand 9 = ir/2 is singled out for the BRT frame, so we chose for a larger value of 8 the value 97r/20. There 
is a symmetry about 9 = 7r/2, so the amplitudes for n/2 < 9 < n don't contain any additional information. 

Eventually we plot the momentum dependence of the valence parts of the form factors and of the violation of the 
angular conditions for two values of the polar angle 9. We see that in all reference frames and for all angles the 
angular condition Abd diverges for Q 2 — » 0. This is due to the fact that the coefficient d,2 of F-i in the matrix element 
G d (See Eq. (|).) vanishes for Q 2 = 0. For that reason there is a finite contribution of the nonvalence part or zero 
mode to the charge form factor even in the limit Q 2 — > 0, which shows up in all d variants of the physical form factors. 
For finite values of Q 2 an accidental singularity in Abd may occur. To follow that up we plot the angular variation 
of the angular conditions for two values of the momentum transfer Q 2 . We can also explain the occurrence of these 
singularities as due to the vanishing of the coefficient di of F-i in the matrix element Gd- 



1. Drell-Yan-West kinematics 



In the DYW reference frame there is no dependence on 9. The dependence on <fi amounts to simple phase factors, 
e ±*4> f or Q b aiK } q c ^ anc j e 2i</> f or Q e j n Figs. ||and||, the results for the valence parts of the invariant form factors 
Fi, F2, and F3 and the physical ones Gc, Gm and Gq are shown. -Fi(O) is normalized to 1, which is not affected by 
the zero mode. For the same reason G b c (0) = G^{0) = 1. However, G^(0) ^ 1 and for positive Q 2 Gq deviates from 
the correct one by the zero-mode contribution to F%. It is clear that the zero mode is very important if one does 
choose the F% -prescription. Neither of the F\- or i^-prescription contains the zero mode. As mentioned before, they 
correspond to the GK-prescription in the DYW-framc for = 0. 



2. Breit frame kinematics 



Our convention for the BRT frame entails both 9- and (^-dependences of the matrix elements. The latter being 
trivial, we fixed <j) — in all our calculations, after checking that indeed the form factors are independent of this angle. 
For 9 = ir/2 the BRT frame and the DYW frame can be connected by a kinematical transformation, so the results for 
the form factors become identical. (See the discussion in ||.) We chose two values for the angle 9, slightly different 
from and tt/2 to illustrate the angular dependence of the valence parts of the form factors. The results shown in 
Figs.|7j-|lC| are for 9 = n/20 and 9 — 97r/20. As shown in Fig.^, it is immediately clear that only F\ and F3 are rather 
insensitive to the choice of the polar angle, but the three prescriptions for F2 are dramatically changing with 9 going 
from a small value to one near 7r/2. This strong angle-dependence is found also in the physical form factors shown 
in Figs.|^-[l0| although the b- and c-variants are much less changed than the d-variant. In all cases the charge form 
factor shows the least angular variation. 
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3. Target-rest-frame kinematics 



The results shown in Figs.|ll|-|14| are again for 9 — n/20 and 9 = 97r/20. Everything we said for the results obtained 
in the Breit frame can be repeated for the target-rest frame. In Figs.[l2|-[l4[ we see for the b- and c-variants similar 
angular dependences, but for the d-variant the variation with 9 is even more dramatic than in the Breit frame. The 
results in FigJlJ for 9 — 97r/20 hint at a singular behaviour of the d- variant that is explained by the fact that f or some 
combinations of Q 2 and 9 the coefficient di vanishes. We discuss more details of the situation below in Section [V B 5. 



4- Angular condition 

In the next plots of Figs Jl5|-|T7|, we show the two angular conditions for a somewhat longer Q 2 interval, up till 10 
GeV 2 . In the case of the Breit (Fig. 16) and target-rest (Fig ,^7|) frames we plot the differences A bc — F$ — F£ and 



be 



A bd = F| - F$ for four angles 9 = tt/8, tt/4, 37r/8,7r/2. For the DYW frame (Fig{l|) A bc = 0, so in the plot A 
coincides with the Q 2 -axis. 

In the B reit frame the angular conditions do depend on 9 and we show their behaviour for the same angles as in 



Sec. IV B 2 . The angular condition that was trivially fulfilled in the DYW frame turns out to be only weakly violated 
in the Breit frame. The other one however, A b d is strongly violated for small values of Q 2 . It demonstrates clearly 
the importance of including the nonvalence parts in a calculation of the matrix elements of the current. For large 
values of Q 2 it tends very quickly to zero, corroborating the expectation that in perturbative QCD one may ignore 
largely the nonvalence parts for sufficiently high momentum tranfers. 

Again, the discussion of the behaviour of the angular conditions in the target-rest frame can be very similar to the 
one for the Breit frame, so we shall not repeat it. We only mention that the overall behaviour is similar in these two 



cases, but the details differ. In particular, in the following subsection 1VB 5| we consider two Q 2 values (Q 2 = 1.0 



GeV 2 and 10.0 GeV 2 ) and find that the singularity in Am occurs in the Breit frame for 9 close to tt/4 while in the 
target-rest frame it shows up for 9 close to 3ir/8. 



5. Singular behavior 

Figure |l4| shows that Gf, starts to drop significantly when Q 2 is near to 4 GeV 2 . Such a behavior can be understood 
from the dependence of the coefficient c? 2 occurring in Eqs. (|4|) and (^), on the momentum transfer Q and the angle 
9. It appears that both in the BRT frame and the TRF this coefficient may vanish for a particular combination of 
Q 2 and 9. 

The singularity of A b d in the Breit frame is illustrated in Fig. [l^. We see that it occurs for Q 2 = 1.0 GeV 2 close to 
9 = 7r/4. The other angular condition remains flat in 9 and the same is true for both conditions for Q 2 = 10.0 GeV 2 . 
A similar picture is found in Fig. [n] for the target-rest frame, only the position of the singularity being different. 

We can understand this behaviour very easily if we consider the expressions for d 2 , that can be derived in a 
straightforward way from the kinematics and the expressions for the polarization vectors inserted in Eq. (Q) . We find 

jBrt = V2/3mQ 2 [l- / 3 2 + (l + /3 2 )cos2fl] 

2 4/3 2 ™ 2 - Q 2 cos 2 9 [ ' 



for the Breit frame and 



(2m + K + PQ cos 9){k 2 + 2f3nQ cos 9 + (3 2 Q 2 cos 29) 



2 2^2m{m + k + f3Qcos9) v ; 

for the target-rest frame. 

Solving the equation d 2 = for Q 2 = 1.0 GeV 2 , we find for BRT 9 = 0.698 = 0.222tt and for TRF 9 = 1.181 = 
0.3767T. These angles coincide with the positions of singularities shown in Figsjl^ and [19[ Also, we note that for 9 = ^ 
in the BRT frame there exists a singularity at Q 2 — 8m 2 /(\/2— 1) ~ 11.45 GeV 2 which is not shown in Fig.fL8| due to 
the restricted interval of Q 2 only up to 10 GeV 2 . However, except the tiny region near this singularity position, the 
angular condition is very well satisfied at the larger Q 2 region. 
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6. Charge, Magnetic and Quadrupole Radii 



From the slope of physical form factors for Q 2 — ► 0, the corresponding radii ({tq), {r 2 M ), (rg)) can be defined^ as 

1 



G C (Q 2 ) ~ G c (0) 
G M (Q 2 ) ~ Gm(0) 
Gq(Q 2 ) ~ Gq(0) 



l-- 6 (r 2 c)Q 2 



1 



1 

To 



r 2 M )Q' 



(44) 



When one considers only the valence parts of charge, magnetic and quadrupole form factors, one should be careful 
in obtaining the corresponding radii given by Eq.(44) because some of them exhibit singular behaviors as Q 2 

G(0))/Q 2 as Q 2 



In order to determine the radii, one may in genera 



0. 

In practice, 



try to take the limit (G(Q 2 ) - G(0))/Q 2 as Q z -► 
however, this gives a rather unreliable value as for very small values of Q 2 the calculations may have numerical noise 
that is amplified by taking the difference of two almost equal numbers and dividing the result by the small number 
Q 2 . A more stable procedure is to make a linear fit to the form factors in a domain close to Q 2 = 0. We have chosen 
the interval 0.01 < Q 2 < 0.1 and took ten equidistant values for Q 2 . In order to check whether it makes sense to fit 
the form factors to a linear function of Q 2 in this domain, we plotted the form factors as shown in Figsj20|-p6| and also 
checked the quality of the fit. 

It turned out that only Gm and Gq in the d-variants could not be fitted with a straight line. The reason is that 
in these variants the influence of the zero mode is very big and only including the zero mode or the nonvalence part 
can correct the nonlinear behaviour. Because of this reason, we do not show the ci-variant valence results for Gm 
and Gq which are anyway out of scale in Figs j2(i|,|23| a nd p6[ However, in the case that includes the nonvalence part 
as shown in Figsj2l|-p3| (total in BRT) and Figs. p4[j2q (total in TRF), all values obtained for the radii do agree. The 
same is true for the case where the zero mode does not occur as shown in Fig.^ (b- and c- variants in DYW). An 
indication of the accuracy of the results is obtained if one includes the values at Q 2 — in the fit. Then the radii 
do not change by more than 1.5%. If one would try do determine these quantities by fitting the form factors in a 
much smaller interval, say .001 < Q 2 < 0.01 in order to improve the linear approximation mathematically, then the 
numerical noise will have a stronger influence. So there is a trade off between truncation error in the series expansion 
of the Bessel function and numerical noise. We are satisfied with an overall numerical error of the order of 1%. In 
Table |, we summarize the numerical results of radii. For the d-variants, both Gm and Gq diverges for Q 2 — > as 
discussed above and an entry 'div.' is given for those cases in this Table. The numerical estimates for the radii in 
these divergent cases, which do not include Q 2 — 0, give indeed values of the order of a hundred GeV -2 . 

TABLE I. Squared radii in GeV -2 for the different reference frames and variants. 



Ref. frame, variant 


Gc 


Gm 


Gq 


DYW, G" and G c 


7.63 


9.73 




12.6 


DYW, G d 


9.56 


div. 




div. 




6 = tv/20 


BRT, tot 


7.63 


9.73 




12.6 


BRT, G b , val 


14.3 


27.6 




50.4 


BRT, G c , val 


14.7 


18.3 




19.4 


BRT, G d , val 


13.9 


div. 




div. 




6 = 97r/20 


BRT, tot 


7.63 


9.73 




12.6 


BRT, G b , val 


8.42 


12.0 




17.7 


BRT, G c , val 


8.48 


10.7 




13.6 


BRT, G d , val 


10.0 
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TRF, tot 


7.63 


9.73 




12.6 



These relations are associated with the interpretation of the Fourier transform of the form factors for spacelike momentum 



transfer as densities and the behaviour of the spherical Bessel functions for small argument ji (x) 
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V. CONCLUSION 



In this work, we made a taxonomical analysis of spin-1 form factors with respect to several different prescriptions 
(GK,CCKP,BH), polarization vector choices (LF,IF), and reference frames (DYW,BRT,TRF) . We used the J+ current 
for all of our analysis. 

In the q + — frame, we looked at both LF and IF polarization vectors and made a comparative analysis on the 
three prescriptions (GK,CCKP,BH) in relating the matrix elements to the physical form factors. In the light-front 
gauge, A + = 0, the light-front helicity basis is the set of eigenvectors of the light-front helicity operator. However, 
the instant-form polarization vectors have been also used in the literature. We find that the zero-mode contamination 
occurs minimally in the light-front gauge because only the helicity zero to zero amplitude, i.e. Gj , gets a zero- mode 
contribution, but all others (G++, Gi , G~+_) are immune from the zero- mode in the light-front basis. Thus, one 
can find a prescription which doesn't involve the zero mode contribution at all. Indeed, the GK prescription which 
doesn't use Gq Q has precisely this property. Consequently, the computation of only the valence contributions to the 
form factors using the GKprescription yields results which coincide exactly with the full results of the form factors 
as we have shown in Figs.gjj. These full results can be obtained by the other prescriptions, CCKP and BH, only if 
the zero-mode contributions are added to the valence contributions. We have also computed the form factors using 
a manifestly covariant Feynman method and explicitly shown that the full results of the light-front calculations are 
fully in agreement with the covariant one no matter what prescriptions we use. Although the full results are also 
independent of the choice of the polarization vectors, we find that the nice feature of the GK prescription described 
above is lost in the instant-form basis. As shown in Figj3[ the valence contribution to the magnetic form factor 
(Gm(Q 2 )) computed in the instant-form basis substantially differs from the full result even if the GK prescription 
is used. This is due to the fact that a larger number of matrix elements are contaminated by the zero-mode in the 
instant -form basis and in the particular case of Gn/fQ 2 ) the terms affected by the zero-mode such as — G+ z and G+,/^/77 
in Eq. (Al) do not cancel each other (See Eq. ( |A6| ).) but rather add up. We thus conclude that the GK prescription 
in the light-front basis is certainly advantageous for model calculation involving only the valence contributions to the 
spin-1 form factors. 

We have also analyzed the frame dependence of the valence contribution to the physical form factors and the 
angular conditions using the light-front polarization vectors. Since the three presciptions discussed above are defined 
only in the q + = frame, we use the prescriptions b, c, and d defined in Section [il] to work in a general frame. In 
DYW frame, our c prescription corresponds to the GK prescription and the results on the form factors from our b 
prescription coincide with those from the c or GK prescription. Also, some combinations of the c and d prescriptions 
correspond to the CCKP and BH prescriptions depending on the coefficients of the combinations. Again, only the d 
prescription involves the zero-mode and the valence result from the d prescription differs from the results of the b and 
c prescriptions that coincide each other exactly as shown in Figs. |5| and |(| The results in the Breit frame reproduce 
the DYW results if 9 — tt/2, since they can be transformed into each other by purely kinematic operators in LFD. If 
9 =/= 7r/2, however, the results are quite different from the DYW results. As shown in Figs. |7[jlC| more drastic differences 
in the results among the 6, c, and d prescriptions are found at 9 = 7r/20 than at 9 = 97r/20. Similar observations can 
be made also for the TRF. However, the kinematic equivalence to DWY obtains only at a special angle 9 = 9q (See 
Ref. |l5|.) which depends on Q 2 . Thus, it is rather difficult to see the similarity of the results of TRF and DYW. 
Although the angular condition A{, c is rather well satisfied, the usual angular condition is severely broken in the 
small Q 2 region. Singularities associated with those violations are visible in our figures(See e.g. Figs.|l8| and p^[). 

Nevertheless, both angular conditions Af, c and A^ are well satisfie d in the region Q 2 above a few GeV 2 except the 
tiny region near the singularity position discussed in Section IV B 5 . and thus the results from all the prescriptions 
become consistent with each other. Thus, one may conclude that the zero-mode contributions are highly suppressed 
in the high Q 2 region and the results are consistent with the perturbative QCD predictions. This may justify the use 
of the BH prescription in the analysis of high Q 2 form factors. For the low and intermediate Q 2 regions, however, the 
zero-mode contributions are very important and the GK prescription with the light-front polarization vectors in the 
DYW frame is certainly desirable for the form factor analyses. An application of this observation to a more realistic 
model calculation is under consideration. 
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APPENDIX A: FORM FACTORS IN THE INSTANT-FORM BASIS 



Using the instant-form linear polarization vectors (h, hi — x, y, z), the form factors corresponding to Eqs. (|l^)-(|l9|) 
in the light-front polarization vectors are given by jq] 
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and 
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(A3) 



where we redefine the definition of the form factors Go, Gi, and G2 in Ref. || in terms of Gc, Gm, and Gq according 
to footnote 1. In order to calculate the form factors of a spin-1 particle in this instant-form bases, the authors 
of Ref. |t],|| used the reference frame specified in our Eq. (|l3|). However, they use different (i.e. non LF gauge) 
polarization vectors in the initial(e) and final(e') states that are given by 



{-VV, Vt + V> 0: 0), e£ = (0, 0, 1, 0), e£ = (0, 0, 0, 1), 



and 



yi^z c z ■ 



(A4) 



(A5) 



(Here we use the component convention p M = (p° , p 1 , p 2 , p 3 ) etc.) Even though these polarization vectors satisfy the 
correct orthonormality and closure relations as well as the conditions e • p = e' • p' = 0, they cannot be obtained in 
the LF gauge. 

Proceeding to calculate the trace terms S£ x , 5+ , S+ and S+. using Eqs. (A4), and (A5), we obtain 



to 2 . + k\ — rjm z 



+ 4k p (p — k )(1 — x — r\ — r\x) 
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4(x + r\x - 1)A£ - 7/(1 + r))x 6 m 2 - (x + 3rjx - 2){m 2 + kj_) 



Q 

yy 

<?+ 



-Ak-{k + -p+)*-4p 4 



4(1 - x)k 2 + r/xm 2 - (2 - x){m 2 + fejj 



= 4fc" 



m„ + fcj_ — r \ m ' L 



+ 4(1 - x)^" 1 



x(l - x)(l + i])m 2 + 2(m 2 + k\) 



-4k 



2 i j 2 2 

m„ + Kj_ — rym 



+ 4y/rjk p + (k + - p + ) 



2k 2 x -x 2 {\- x){\ + r])m 2 - 2{m 2 + k\) 



(A6) 



where ^2p + = p° = rriy/1 + rj and we omit the terms of odd power in k± since they do not contribute to the integral. 
Note also that we do not separate the on-shell propagating part from the instantaneous one in this instant-form 
calculation. As we discussed, only the underlined terms in S xx , Sf z and S+ x contribute to the zero- mode part and the 
S+ y component is immune to the zero-mode contributions. 



APPENDIX B: DERIVATION OF THE ZERO-MODE 



In this appendix, we derive Eq. ( |39| ) in detail. Performing the k -integration in Eq. (|39|), the residue at the pole 
k- = p'+ + [(k ± - p\f + A 2 }/2(k+ - p'+) gives 
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(-q+)(k+ -p+)(k+ -p'+)(m 2 +p\) +p+p'+(k+ -p+)[(k± -p' ± ) 2 + A 2 ] 
-p + p' + (k + - p' + )[(kj_ - p±) 2 + A 2 ] 



(Bl) 



where we use p]_ = p'j_ and q + = p' + — p + in the derivation of the last term in Eq. (Bl). Writing q + = 5p + and 
p l+ = (1 + S)p + , 5 — > at the end, we obtain 
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If we write x = 1 + S y and dx — S dy, then the integral over y runs from to 1 as x runs from 1 to 1 + S. Therefore, 
we get 
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FIG. 2. The charge form factor \Gc{Q 2 )\ obtained from the light-front (left) and the instant-form (right) spin bases: The 
thick solid line represents the full (i.e. valence+zero-mode in LF = covariant) solution. The dotted, long-dashed, and dot-dashed 
lines represent the valence contributions only, where we use the same normalization as for the full solution Gc(0) = 1. The 
small squares represent the angular condition in Eq. (|l^) without including the zero-mode contribution. 
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FIG. 3. The magnetic form factor Gm(Q 2 ) obtained from the light-front (left) and the instant form (right) spin bases. The 
same lines are used as in Fig 0. 
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FIG. 4. The quadrupole form factor Gq(Q 2 ) obtained from the light-front (left) and the instant form (right) spin bases. The 
same lines are used as in Fig H. 
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FIG. 5. Invariant form factors Fi, F$, F%, F$, and F3 (left) and physical form factors Gc, Gm, and Gq from F 2 (right) 
calculated in the Drell-Yan-West frame. Valence parts only. 
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FIG. 9. Physical form factors Gc, Gm, and Gq from F 2 C calculated in the Breit frame. Left 9 = it/20, right 9 = 97r/20. 
Valence parts only. 
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FIG. 10. Physical form factors Gc, Gm, and Gq from F$ calculated in the Breit frame. Left 8 — tt/20, right 8 — 97r/20. 
Valence parts only. 
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FIG. 12. Physical form factors Gc, Gm, and Gq from F\ calculated in the target-rest-frame. Version b. Left 8 = 7r/20, 
right 9 = 97r/20. Valence parts only. 
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FIG. 13. Physical form factors Gc, Gm, and Gq from F 2 C calculated in the target-rest-frame. Left 6 = n/20, right # = 97r/20. 
Valence parts only. 
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FIG. 16. Dependence of the angular conditions Ab c and Abd in the Breit frame on Q for four different angles 9. 
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FIG. 17. Dependence of the angular conditions Ab c and Abd in the target-rest frame on Q 2 for four different angles 9. Note 
the changed scale in the case 9 = 37r/8. The singlularity is clearly visible there. 
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FIG. 18. Angular conditions Abe and Abd in the Breit frame for Q 2 = 1.0 GeV 2 and Q 2 = 10.0 GeV 2 for different angles 
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FIG. 19. Angular conditions Abe and Abd in the target-rest frame for Q 2 = 1.0 GeV 2 and Q 2 — 10.0 GeV 2 for different 
angles 6. 
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FIG. 20. Valence contributions to physical form factors for small values of Q in DYW frame. While b and d variants are 
shown here, the d-variant magnetic(G^) and quadrupole(GQ) form factors are out of scale because they diverge as Q 2 — > 0. 
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FIG. 22. Physical form factors for small values of Q 2 . BRT frame, variant c. 
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FIG. 23. Physical form factors for small values of Q 2 . BRT frame, variant d. 
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FIG. 24. Physical form factors for small values of Q 2 . TRF, variant b 
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FIG. 25. Physical form factors for small values of Q 2 . TRF, variant c. 
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FIG. 26. Physical form factors for small values of Q 2 . TRF, variant d. 
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